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Abstract. We discuss our ongoing program of Target of 
Opportunity sub-millimeter observations of gamma-ray 
bursts (GRBs) using the Sub-millimetre Common-User 
Bolometer Array (SCUBA) on the James Clerk Maxwell 
Telescope (JCMT). Sub-millimeter observations of the 
early afterglows are of interest because this is where the 
emission peaks in some bursts in the days to weeks follow- 
ing the burst. Of increasing interest is to look for under- 
lying quiescent sub-millimeter sources that may be dusty 
■ star-forming host galaxies. In this paper, we present obser- 
vations of GRB 991208, 991216, 000301C, 000630, 000911, 
and 000926. For all these bursts, any sub-millimeter emis- 
sion is consistent with coming from the afterglow. This 
means that we did not conclusively detect quiescent sub- 
millimeter counterparts to any of the bursts that were 
studied from 1997 through 2000. The inferred star forma- 
tion rates (M > 5Mq) are typically ^ 300 M0yr-^ If 
GRBs are due to the explosions of high-mass stars, this 
may indicate that the relatively small population of ex- 
tremely luminous dusty galaxies does not dominate the 
total star formation in the universe at early epochs. In- 
stead, the GRBs may be predominantly tracing slightly 
lower luminosity galaxies. The optical faintness of some 
host galaxies is unlikely to be explained as due to dust 
absorption in the host. 

Key words: gamma rays: bursts - radio continuum: gen- 
eral - infrared: general 
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1. Introduction 

The discovery of localized transients in the error boxes of 
gamma-ray burst (GRB) sources has led to intense multi- 
wavelength campaigns that have revolutionized our under- 
standing of these so urces . For reviews, see Piran ( p^99| ); 
Van Paradijs et al. (^000|) . 

The multi-wavelength emission comes from several 
components. The "prompt" emission comes from the ini- 
tial explosion. A reverse shock can give an optical and/or 
radio flash. The later "afterglow" emission comes from 
the expanding fireball as it sweeps up the surrounding 
medium. Finally, the "quiescent" constant emission comes 
from any underlying host galaxy. 

Sub-millimeter observations provide "clean" measures 
of the source intensity, unaffected by scintillation and ex- 
tinction. We have been performing Target of Opportu- 
nity sub-millimeter observations of GRB counterparts us- 
ing the Sub-millimetre Common-User Bolometer Array 
(SCUBA) on the James Clerk Maxwell Telescope (JCMT) 
on Mauna Kea, Hawaii. The detailed SCUBA results for 
the first eight (GRB 970508, 971214, 980326, 980329, 
980519, 9807 03, 9 81220, and 981226) are described in 
Smith et al ( |1999D . GRB 990123 is discussed in Galama 
et al. ( |1999| ). Observations of G RB 99 0520 were made in 
mediocre weather (Smith et al. 2000a). In this paper we 



present our SCUBA observations of GRB 991208, 991216, 
000301C, 000630, 000911, and 000926. The sub-millimeter 
emission from GRB 010222 that was detected b y SCU BA 
appeared to ha ve unusual properties (Fich et al. 2001 ; Ivi- 
son et al. 2001 ); these observations will be reported else- 
where after a final sensitive SCUBA observation has been 



made to determine whether there is an underlying host 
galaxy. 

In the remainder of §1, we outline the motivations for 
making sub-millimeter observations of the afterglows and 
host galaxies. In §2 we discuss the technical details of the 
SCUBA observations, and also the chance probability of 
detecting quiescent sources that are unrelated to the GRB. 
In §3 we present the results of our latest observations. In 
§4 we discuss some of the implications of our results. 

1.1. SCUBA afterglow ■program 



Bo th ob servations (e.g. Galama et al. 1998 , 199£ 
al 



Smi th et 

and th eories (e.g. Sari et a l |1998| ; Piran |l999| ; Wi- 
jers & Galama 1999; Granot et al. 200C; Sari & Meszaros 
^OOOt Chevalier & Li |2000| ; Panaitescu & Kumar ^OOOj ) 
show that the burst afterglow emission often peaks in the 
sub- millimeter in the days to weeks following the burst. 
By tracking the evolving emission across the entire spec- 
trum, it is possible to study aspects such as the types of 
shocks involved, the geometry of the outflow (jet versus 
spherical), and the geometry of the surrounding medium 
(uniform versus stellar wind). 

Previously, our most interesting result was the detec- 
tion of a fading sub-millimeter counterpart to GRB 980329 
(Smith et al. 1999). While a fading X-ray counterpart (in't 
Zand et al. 1998) and a variable, long-lasting radio source 
VLA J070238.0+385044 (Taylor et al. |1998D were soon 



found for this burst, it proved to be difficult to find the 



infrared counterpart (Klose et al. 1998; Palazzi et al. 1998 
Reichart et al. |1999| ). However, the 850 ^im (350 GHz) flux 
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seen by SCUBA was relatively bright, making it similar 
to GRB 970508. The radio through sub-millimeter spec- 
trum of GRB 980329 was well fit by a power law with a 
surprisingly steep index of a = -1-0.9 (Taylor et al. 1998; 
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Smith et al. |1999D . However, we could not fully exclude a 
v^l^ power law attenuated by synchrotron self-absorption 



(e.g. Katz |l994t Waxmau pM?! ). The 850 /^m SCUBA flux 
decayed rapidly with time. For a power law decay with the 
flux density oc t~™ where t was the time since the burst, 
the best fit power law index was m = 3.0. However, m was 
not tightly constrained: the 90% confidence interval was 
m — 1.2 to m = 5.3. 

1.2. SCUBA host galaxy program 

It is important to understand the nature of the host galax- 
ies. A popular model for GRBs is that they are due to the 



explosive deaths of high-mass stars (e.g. Woosley 1993 



200C; Paczyhski 1998) and therefore are likely to be found 



in active star forming regions. SCUBA has recently dis- 
covered several dusty sta r-forming g alaxies out to high 
redshifts (e. g. Sma il et al. |l997| , |l998l Hughes et al. |l998 ; 
Barger et al. 1998 ), and it appears that the star formation 
rate does not drop rapidly beyond z ~ 1 (e.g. Blain et al. 
1999). In this scenario, it is plausible that some of the 



GRBs will be associated with dusty star- forming galaxies. 

An alternative model assumes that GRBs are due to 
the coalescence of two compact objects that were in a bi- 
nary orbit. In this case, the lifetime of the binary will likely 
be long, and the active star formation will have died down 
by the time the burst occurs. Furthermore, the binary sys- 
tem may have obtained a substantial velocity if there was 
a supernova explosion with a large kick; this would move 
the system far from its birthplace. Thus GRBs would be 
less likely to occur in a dusty star-forming galaxy. 

If there is a connection between GRBs and dust- 



enshrouded star formation (Blain & Natarajan 200C), it 
has been suggested that ^ 20% of GRB hosts sho uld b e 
brighter than 2 mJy at 850 /xm (Ramirez- Ruiz et al. 2001). 
Our Target of Opportunity program is designed to look for 
the afterglow emission of the bursts by making observa- 
tions separated by days. However, by combining the data 
from all of our observations of a source we can also look 
for quiescent sources of this nature. 

The expected sub-millimeter emission from warm dust 
due to massive star formation is somewhat uncertain. In 
this p aper, we will gi ve nu merical estimates using (Condon 
199^ ; Carilh & Yun [19991) : 



5^ = 0.039 ^(z,go,a) 



Ha 



65 kms Mpc 



L350 GHz J 



SFR (M>5Mq) 
100 MQyr-i 



mJy , 



(1) 



where Sy is the sub-millimctcr flux density at the fre- 
quency v. 

SFR is the star formation rate for stars with masses 
> 5Mq in units of 100 MQyr"-'^. An advantage of using 
the sub-millimeter flux density to estimate the SFR is that 
there is no need to perform the large (sometimes orders 
of magnitude) correction for the host galaxy extinction 
that is required when studying optical lines. Note that 



other methods of estimating the SFR often quote results 
using masses between 0.1 and 100 Mq. Our results will 
be smaller than these by a factor that depends on the 
initial mass function (IMF). This is currently uncertain 
for starburst galaxies; for a Salpeter IMF, the factor is 
- 5. 

The function ^(z, q^, a) is given by: 
(H-z)i+" 



f(z,go,a) = 



1 I ^(i-go) 

^ (l+2goz)i/^ + l+go2 



(2) 



This assumes that the Universe has a geometry described 
by the Robertson- Walker metric and uses the Friedmann 
cosmological models. Typically, i,{z,qi),a) ^ 10 over a 
wide range of redshifts, assuming plausible values for go 
and a. In our numerical calculations, we assume Hq = 
65 kms~^Mpc"\ flo = 0.2, qo = 0.1, and Aq = 0. 

The power law index a is measured between 230 and 
850 GHz, and is expected to be 3 — 4 based on observa- 
tions and simple dust models (e.g. Dwek & Werner 1981). 
Thus any quiescent dust contribution is very much larger 
at sub-millimeter than at radio wavelengths. In principle, 
a could be measured by SCUBA. However, in practice 
the weather is rarely good enough to be able to detect the 
source at 450 /im. For our GRB numerical calculations, 
we assume a = 3.4 (which is the value observed for M82). 
This lies in the middle of the expected range of a, and 
for the z < 2 bursts considered here the uncertainty in a 
means that the calculated SFR may be too high or too 
low by at most 50%. Note that depending on the tem- 
perature of the dust, for very large z > 7 the peak of the 
dust emitting spectrum may be redshifted to a wavelength 
^ 850 fim; in this case, a will be much flatter. 

As an example, a starburst galaxy with a SFR = 
300 MQyi~-^, and a = 3.0 (as in Arp 220) would give 
an 850 /im flux density of Si, = 1.5 niJy at z = 10. Al- 
ternatively, a flux density of S^, = 1.5 mJy at 850 /im for 
an ultraluminous infrared galaxy (ULIRG) at z = 1 with 
a = 3.4 would imply a SFR = 360 MQyr^^. 

SCUBA is therefore capable of detecting ULIRG-type 
objects throughout the distant universe. The relatively 
small population of extremely luminous dusty galaxies 
may dominate the total star formation in the universe at 
early epochs (Smail et al. 2000 ; Ramirez- Ruiz et al. 2001). 
However, this is still quite uncertain; for example, there 
may be a substantial contribution to the energy output of 
the SCUBA-bright galaxies from active galactic nuclei. If 
GRBs are regularly found to be associated with extremely 
luminous dusty galaxies, this will provide independent ev- 
idence that these objects dominate the star formation in 
the early universe. 

In Smith et al. (1999) we noted that the relatively 
bright sub-millimeter flux from GRB 980329 could be due 
in part to an underlying quiescent sub-millimeter source. 
Unfortunately, the quiescent flux of 1 mJy at 850 /im 
in this case is too faint to be significantly detected by 
SCUBA: more sensitive sub-millimeter instruments will 
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be needed to determine whether such a quiescent source 
is actuahy present in GRB 980329. 

None of the other GRBs presented in our previous pa- 
pers showed any indication of an underlying quiescent sub- 
miUimeter source. We discuss the SFR limits for these in 
§4. However, 60 fim Infrared Space Observatory observa- 
tions suggested a possible far-infrared quiescent counter- 
part to GRB 970508 (Hanlon et al. ^OOOj). 



2. SCUBA details 

SCUBA is the sub-millimeter continuum instrument for 
the JCMT (for a review see Holland et al. |1999D . Our 
observing, calibration, and reduction techniques are the 
same as described in detail in Smith et al. ( |1999D . Here we 
summarize only the most important features of the instru- 
ment, and those aspects that have been modified since our 
previous observations. We also discuss the chance proba- 
bility of detecting a quiescent source that is unrelated to 
the GRB. 



2.1. SCUBA observing details 

SCUBA uses two arrays of bolometers to simultaneously 
observe the same region of sky, ~ 2.3' in diameter. The 
arrays are optimized for operations at 450 and 850 /im. 
The dedicated photometry pixels for observations at 1100, 
1350, and 2000 /im have not been available since the 1999 
upgrade. 

In 1999 October new wideband filters (450w:850w) 
were installed. These have improved sensitivities over the 
ones used previously by SCUBA. New blocking filters with 
better transmission were also installed in 1999 October. 
New ribbon cables now allow the refrigerator to run colder: 
this has reduced the number of noisy pixels and has sta- 
bilized the noise on the arrays. This updated set-up was 
used for all the new observations described here, and has 
generally resulted in lower rms values than was found in 
Smith et al. ( |1999D . A typical integration time of 2 hours 
gives an rms ~ 1.5 mJy at 850 /im. However, the sensitiv- 
ity depends significantly on the weather and the elevation 
of the source; since our ToO observations are done on short 
notice, sometimes these factors are less than ideal. 

Fully sampled maps of the 2.3' region can be made by 
"jiggling" the array. However, for all the sources described 
here, we have only been looking at the well-localized op- 
tical or radio transient coordinates by performing deep 
photometry observations using a single pixel of the arrays. 
The other bolometers in the arrays are used to perform a 
good sky noise subtraction. 

During an observation the secondary is chopped be- 
tween the source and sky at 7 Hz. This is done mainly 
to take out small relative DC drifts between the bolome- 
ters, and also to remove any large-scale sky variations. 
The term "integration time" always refers to the "on-|-off" 
time, including the amount of time spent while chopped 



off-source. An 18 sec integration thus amounts to a 9 sec 
on-source observation time. A typical measurement con- 
sists of 50 integrations of 18 seconds; we refer to this as a 
"run." Each observation of a source in general consists of 
several such runs, with pointing and calibration observa- 
tions in between. 

Based on observed variations of the gain factor and 
signal levels we estimate typical systematic uncertainties 
in the absolute flux calibrations of 10% at 850 //m. In 
general the rms errors of the observations presented here 
are larger than this uncertainty. 



2.2. False positives 

The pointing accuracy of the JCMT is a few arcsec, and 
the pointing is checked several times during the night to 
ensure that it is reliable. The 850 /xm bolometric pixel has 
a diffraction limited resolution of 14". Thus we can be sure 
that the target is always well centered in the bolometric 
pixel. 

However, the large beam size combined with the large 
number of distant galaxies radiating strongly at this wave- 
length means that in any observation there is a non- 
negligible chance of detecting a quiescent sub-millimeter 
source that is completely unrelated to the GRB. 

The surface density of sub- millim eter galaxies is sti ll 
somewhat uncertain (Blain et al. 199^ ; Barger et al. 1999). 
At 850 fiui, the surface density of galaxies with flux den- 
sities larger than 4 mJy is estimated to be ~ (1 — 2.5) x 
10'^ deg~^, while the surface density of galaxies with flux 
densities larger than 1 mJy is ^ 10^ deg~^. 

Because of the uncertainty in the surface density, and 
the fact that the sensitivity of the bolometric pixel de- 
pends on the off-axis distance of the source, it is difficult 
to give an accurate determination of the chance probabil- 
ity of detecting a quiescent sub-millimeter source that is 
completely unrelated to the GRB. However, we estimate 
that the chance of detecting a random > 4 mJy source in 
any pointing is ~ 1 — 3%, while the chance of detecting a 
random > 1 mJy source in any pointing is > 10%. 

The rather large chance of a false positive means that 
it is dangerous to use just the detection of a quiescent 
sub-millimeter source to claim that this must be the host 
galaxy to the GRB. For example, we cannot rule out that 
the 1 mJy source that may have been seen in GRB 980329 
is an unrelated source in the JCMT beam. For individual 
cases, confirmation of the star formation rate is needed 
from observations at other wavelengths. However, if it is 
found that many more bursts are associated with quiescent 
sub-millimeter sources than is expected by chance, this 
would be good evidence that the majority of these are 
true associations. 
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Table 1. SCUBA 850 fim (353 GHz) GRB observations. 



Jjurstr 


Z 


I 1 i~\c?i^"V"i ri Ti rr 
KJ Ubci V iligi 

date 
(UT) 


j_ line hiiice 

burst 

(days) 


iim fliiY 
OOU ^111 iillA 

density 
(mJy) 


991208 


0.707 


19991215.82 


7.63 


3.4 ±3.7 






19991219.82 


11.63 


-0.8 ± 1.8 


991216 


1.0 


19991218.48 


1.81 


0.7 ± 1.6 






19991219.45 


2.78 


-2.0 ± 1.7 


000301C 


2.04 


20000304.75 


3.34 


3.1 ±3.1 






20000305.53 


4.12 


1.9 ± 1.2^ 






20000306.50 


5.09 


1.1 ±0.9 


000911 




20000917.53 


6.23 


-0.4 ± 1.4 






20000920.49 


9.19 


0.3 ± 1.1 


000926 


2.04 


20000930.27 


3.28 


7.3 ±4.2 



V = 24.6 magnitude galaxy was later found at this loca- 
tion (Fruchter et al. piOOOD . 

Unfortunately, GRB 991208 was in the morning sky 
and during a period of very poor weather at Mauna Kea. 
Our first opportunity to make a (short) observation was 
not until a week after the burst. Table 1 shows the re- 
sults of our two observations. They are consistent with 



' This result uses all of our data, including a negative out- 
lier. Alternatively, we could quote a "best" result of 3.3±1.2 
mJy around the time of the peak magnification in the lens- 
ing scenario. See the text for details. 



the other millimeter results (Galama et al. 200C), though 
we are unable to determine anything significant about the 
evolution of the flux. 

A SFR of (11.5 ± 7.1) M0yr-i was estimated from 
the e missio n line fluxes of the host galaxy (Castro-Tirado 
et al. 2001). The expected quiescent flux density at 850 
fj,m would then be 5^ ~ 0.06 mJy, assuming a = 3.4. 
However, correcting for the internal extinction in the 
host galaxy leads to significantly higher SFR estimates, 
> 156 Mgyr-i to > 249 Mgyr"! (Sokolov et al. ^OOlj ). 
A SFR {M > 5Mq) of 250 Mqyt-^ would give a flux 
density of ~ 1.2 mJy at 850 fim; since we do not detect 
a quiescent sub-millimeter source with SCUBA, the SFR 
cannot be significantly higher than this value. 



3. Results of SCUBA observations 

The fact that SCUBA was down during the 1999 upgrade 
combined with unusually poor weather and a lack of well- 
localized sources in regions of the sky accessible to SCUBA 
has unfortunately limited our program over the past cou- 
ple of years. For completeness, we describe here all of our 
observations from late 1999 through 2000. Table 1 sum- 
marizes the 850 /im results. 

3.1. GRB 991208 

The extremely energetic GRB 991208 was detected 



19991208.19 (Hurley et al. |2000a| ) by the Ulysses, Near- 
Earth Asteroid Rendezvous (NEAR), and Wmd spacecraft 
of the third Interplanetary Network (IPN). 

A radio source was found by the Very Large Array 
(VLA: Frail |1999|; Hurley et al. |2000a|). This source had 



a rising spectral index. A few days after the burst, the 
source was found to have fluxes of a couple of mJy from 
15 to 240 GHz using the Ryle Telescope (Pooley |l999aD , 
Owens Valley Millimeter Array (Shepherd et al 
and IRAM (Bremer et al. [19991) . 



19991) 



An o ptical transient was a lso found (Castro-Tirado et 



al. |1999| ; Stecklum et al. |1999[ ) whose optical flux decayed 
following a power law c>c i^* with a steep slope S = 2.3 ± 
0.07 for the first ~ 5 days, breaking to an even faster 



(5 = 3.2 ± 0.2 (Castro-Tirado et al. |2001|). 

A redshift of z = 0.707 was determined for the after- 



glow (Dodonov et al. 1999; Djorgovski et al. 1999; Castro 



Tirado et al. 2001; Sokolov et al. 2001). There was an in 



dication of a supernova or dust echo emission component 
in the afterglow (Castro-Tirado et al. ^001 ). A compact 



3.2. GRB 991216 

The extremely bright GRB 991216 was detected 
19991216.67 by BATSE on the Gompton Gamma-Ray Ob- 



servatory {GGRO) (Kippen et al. |1999| Kippen |1999|) , and 
by Ulysses, NEAR, and the Satellite per Astronomia a 
Ragg i X (BeppoSAX) (Hurley & Feroci |1999| ; Hurley et al. 
|1999| ). Scans of the GRB error box using the Proportional 



Counter Array (PGA) on the Rossi X-ray Timing Explorer 
(RXTE) discovered a p revio usly unknown fading X-ray 
source (Takeshima et al. 1999|) . This source was accurately 
locat ed by the Ghandra X-Ray Observatory (GXO) (Piro 
et al. |1999| ). 

The X-ray transient location was consistent with that 
of a fading optical source (Uglesich et al. 1999), whose 
optical power law decay of S 1.2 — 1.5 was fla tter than 
the X- ray d ecay of ^ 1.6 (Garnavich et al. 2000a ; Halpern 
et al. 2000). Detections of both redshifted iron features 



(Piro et al. 2000) and optical absorption lines (Vreeswijk 
et al. 1999) indicated that the burst was at 2 = 1.0. An 
irregularly shaped R = 26.9 magnitude galaxy was later 
found at this location (Vreeswijk et al. 2000| ). 

A variable radio afterglow was also found soon after the 
burst (Taylor & Berger 1999| ; Rol et al. 1999D . Surprisingly, 
the radio decay started within 1.5 days of the burst, and 
the radio power law decay index was less s teep than in 
both the optical and X-ray bands (Frail et al. 2000 ). While 
a 1.1±0.25 mJy source at 15 GHz was see n 1.33 days after 
the burst by the Ryle Telescope (Pooley 1999b| ), this was 
not detected in later Ryle observations, nor was it detected 
at 99.9 GHz by Owens Valley (Frail et al. |2000D . 

Unfortunately, as with GRB 991208, the weather at 
Mauna Kea was very poor during this period. Our SCUBA 
results are listed in Table 1. Our upper limits are consis- 
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tent with the lack of a detection by Owens Valley. We find 
no evidence for a quiescent dusty host galaxy for GRB 
991216. Combining all our 850 /xm runs, the rms value 
of 1.2 mJy implies a SFR (M > SM©) < 290 M^yi^'^, 
assuming a — 3.4. 

3.3. GRB 000301 C 

The short or intermediate duration (2 s in the > 25 keV 
energy range) GRB 000301C was detected 20000301.41 by 
the AU Sky Monitor (ASM) on RXTE, b y Ulys ses, and by 
NEAR (Smith et al. |2000bt Jensen et al. |2001| ). 

An optical counterpart was soon found in the burst 
error box (e.g. Fynbo et al. 2000a; Jensen et al. 2001). 
The light curve of this relatively bright optical source was 
well sampled, and it was found that it had a complex be- 
havior, with significant fluctuations superimposed on the 
overall decay (e.g. Masetti et al. 200C; Berger et al. 2000a). 
These variations might be due to refreshing of the shock 
in the fireball (Dai & Lu 2001), or they may be due to 
inhomogeneities in the ambient medium that the fireball 
is expanding into (Berger et al. 2000a). The achromatic 
fluctuations have also been explained as a micro-lensing 
event due to a 0.5 Mq le ns located half way to the burst 
(Garnavich et al. 2000b; but see also Panaitescu 2001): 
the peak magnification of ~ 2 occurred 3.8 days after the 
burst. 



Detections of a Lyman alpha break (Smette et al. 2001) 
and weak optical absorption lines (Castro et al. 2000a; 
Feng et al. 2000; Jensen et al. [2001 ) indicated that the 
burst was a,t z = 2.04. A faint, extended, i? = 28 mag- 
nitude object was tentatively suggested to be the host 
galaxy (Fruchter & Vreeswijk 2001). 

At the location of the optical transient, a radio and 
millimeter source was found by the VLA, Ryle, Owens 
Valley, and IRAM (Berger ct al. |2000a| ). Our SCUBA re- 
sults are listed in Table 1. Unfortunately, our observations 
on the first day had to be terminated early, resulting in a 
large rms. 

The SCUBA results for the second day are somewhat 
confusing. While it appeared that a source was present, 
one of our six runs had a rather negative result. There is 
some indication that the atmosphere was a little unstable: 
there were a couple of deviations in the opacity monitored 
by the Caltech Sub-millimeter Observatory. However, the 
JCMT was probably looking in a completely different di- 
rection at the time. In Table 1, we have conservatively 
chosen to show the result that uses all of our data, includ- 
ing the negative outlier. If we drop the negative run and 
combine with the data from the first day, we instead get 
a "best" result of 3.3 ±1.2 mJy around the time of the 
peak magnification in the lensing scenario. This is consis- 
tent with the other multi-wavelength results that showed 
that the broad-band spectrum was peakin g in the sub- 
millimeter at this time (Berger et al 



Our long SCUBA observations on the third day gave a 
lower rms, but did not result in a detection of the source at 
the level suggested from the previous day. This is fully con- 
sistent with the results at other wavelengths, that showed 
that the achromatic fluctuation had subsided. However, 
given the weak significance of our results, we caution 
against over-interpreting the possible transient magnifi- 
cation of the 850 /im flux. 

Combining the data from all three days gives a flux 
density of 1.9 ± 0.7 mJy. The rms is the lowest value we 
have obtained for any burst to date. In this case, the flux is 
consistent with coming entirely from the afterglow. Thus 
there is no evidence for a quiescent dusty host galaxy for 
GRB 000301C. Using the 850 fiin. rms value implies a SFR 
(M > 5Mq) < 175 M0yr-\ assuming a = 3.4. Thus 
the optical faintness of the host galaxy is unlikely to be 
explained as due to dust absorption in the host; this is 
consistent with the low extinction derived from the optical 
properties of the afterglow. 

3.4. GRB 000630 

GRB 000630 was detected 20000630. 02 by Ulysses, Wind, 
NEAR, and BeppoSAX (Hurley et al. |2000bD . An unknown 
radio so urce was found inside the GRB error box (Berger 
& Frail 2000). However, a faint fading optical source was 
also f ound in the error b ox at a different location (Jensen 
et al. 2000; Fynbo et al. 2001a ). A faint quiesce nt obj ect 
was found at this optical location (Kaplan et al. 2001). 

Only one SCUBA observation was attempted on this 
burst, on July 2, with the radio source as the target. The 
weather was very poor, and no source was detected with 
an rms of 2 mJy at 850 /zm. Since the optical source was 
the more likely counterpart, we have not included this 
result in Table 1. 



3.5. GRB 000911 

The long GRB 000911 was detected 2 00009 11.30 by 
Ulysses, NEAR, and Wind (Hurley et al. ^000c|) . A new 
radio and optical source was found inside the GRB error 
box (Berger et al. |2000bD . 

The optical power law decay was i5 ^ 1.5 (Price et 
200C| ). In the simple adiabatic pis- 



al. 2000; Lazzati et al. 



ton model, a flreball produced by a one time impulsive 
injection of energy in which only the forward blast wave 
efSciently accelerates particles predicts a power law spec- 
trum Su oc w ith energy spectral index 13 = 25/3 



(Wijers et al. |l997D . For GRB 000911, this would imply 
/3 ^ 1.0. While the measured broad-band optic al spec- 
trum had a steeper /3 ^ 1.5 (Lazzati et al. 2000 Covino 
et al. 2000 ), correcting for Galactic extinction brings the 



2000a) 



spectrum closer to that of the simple fireball model. 

The SCUBA results from our first two observations on 
the radio-optical counterpart are listed in Tabl e 1. T hese 
were originally reported in Smith & Tilanus (2000). We 
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did not detect the source at 850 /zm. A thir d SC UBA 
observation that was listed in Smith & Tilanus ( |2000| ) was 
made on the incorrect coordinates given in Price et al. 



Table 2. Inferred SFR hmits (M > 5Mq) using SCUBA 
850 /im rms values, and a = 3.4. 



(2000: see Price 2000), and is not included here. 

It is currently difficult to meaningfully compare our 
SCUBA limits with the results of the afterglow emission 
at other wavelengths. Simply scaling from the radio data 
using our fit for GRB 980329 would have led to a higher 
850 //m flux than we observed. However, the GRB 000911 
radio data still needs to be corrected for possible interstel- 
lar scintillation effects. Similarly, extrapolating the optical 
spectrum would also have predicted a higher 850 /xm flux 
than we observed, although the final extinction correc- 
tions need to be incorporated. It does appear, however, 
that the SCUBA results will ultimately imply there are at 
least two breaks in the spectrum between 10'^"'^ and 10^^ 
Hz, as is seen in other afterglows. 

Combining all the SCUBA data gives an rms of 0.9 
mjy. Thus there is no evidence for a quiescent dusty host 
for GRB 000911. There has also not yet been any indi- 
cation of a host galaxy in the optical observations. Since 
no redshift has been reported so far, we are unable to put 
an accurate limit on the SFR, although it is likely to be 
< 200 MQyr^i. As for GRB 000301C, the optical faint- 
ness of the host galaxy is unlikely to be explained as due 
to dust absorption in the host. 



3.6. GRB 000926 

GRB 000926 was detected 20000926.99 by Ulysses, Wind 



and NEAR (Hurley et al. |2000d|). BeppoSAXand CXO ob- 
servations located a fading X-ray counterpart in the burst 
error box that had a very steep power law decay (Piro 
2000| ; Pho & Antonelh |2O(30| ; Garmire et al. |2000D . 

Consistent with the location of this X-ray source, 
a bright new optical and radio counterpart was found 
(Gor osabel et al. ^OOOt D aU et al. |2000| ; Frail & Berger 
2000| ; Fynbo et al. 2001b| ). The optical decay showed a 
sharp achromatic steepening with time (e.g. Fynbo et al. 
2000b| , ^001b| ; Rol et al. ^000| ; Price et al. |00l|). An i? - 25 
magnitude compact knot of emission was found at this 
location, which may be the host galaxy (Harrison et al. 
2001). A redshift of z = 2.04 was determined using the 



optical absorption lines in the afterglow spectrum (Fynbo 



et al. 2000c; Castro et al. 2000b) 



Unfortunately, GRB 000926 was in the afternoon sky 
during a period of poor weather at Mauna Kea. Only one 
long SCUBA observation was performed, in very marginal 
weather. The result is listed in Table 1. While the result 
was tantalizing, it was not possible to obtain a confirma- 
tion. Extrapolating from the 0.5 mJy flux density at 
8.46 GHz (Frail & Berger |200C| ) , the 850 /im flux density 
would be expected to be 1.7 mJy using a 5',^ oc v^^^ 
spectrum, or - 14 mJy for oc i'^-^ (as in GRB 980329; 
see §1.1). Thus the SCUBA flux is consistent with coming 
from the afterglow. 



Burst 


z 


850 fim rms 
(mJy) 


SFR (M > 5Mq) 
(M0yr-i) 


971214 


3.418 


1.0 


< 200 


980519 




1.8 




980703 


0.966 


1.6 


< 380 


990123 


1.6 


0.7 


< 180 


991208 


0.707 


1.8 


< 370 


991216 


1.0 


1.2 


< 290 


00030 IC 


2.04 


0.7 


< 175 


000911 




0.9 





The SFR was estimated to be 14 MQyr"^ without 
including extinction in the h ost ga laxy, and 24 MQyr~^ 
with extinction (Fynbo et al. 2001b). Using the extinction 
corrected value, the expected quiescent flux density at 850 
/im would then be < 0.1 mJy, assuming a = 3.4. Thus we 
would not expect to significantly detect a quiescent sub- 
millimeter source using SCUBA. 



4. Discussion 

With the possible exception of GRB 980329, for aU of 
the bursts that were observed by SCUBA between 1997 
and 2000, any sub- millimeter emission is consistent with 
coming from the afterglow: we did not conclusively de- 
tect quiescent sub-millimeter counterparts for any of these 
sources. 

Unfortunately, the redshift remains uncertain for GRB 
980329. If we assume that there is a 1 mJy 850 /xm quies- 
cent source present, then we obtain a SFR (M > 5Mq) 
of 240 Mqyt-^ if z = 1 (using a = 3.4), or 160 Mqjt-'^ 
if z = 5 (Fruchter 1999] ). However, since the chance of 
detecting a random > 1 mJy source in any pointing is 
> 10%, we cannot rule out that the 1 mJy quiescent emis- 
sion that may have been seen in GRB 980329 is from an 
unrelated source in the JCMT beam. 

Table 2 summarizes the SFR limits for the other bursts 
that have had secure counterpart determinations and that 
had rms values less than 2 mJy measured by SCUBA. For 
those bursts with reasonably reliable redshift determina- 
tions, the SFR limit has been estimated. Our sample spans 
a wide range of redshifts. No redshifts have so far been re- 
ported for GRB 980519 or GRB 000911; however, since 
the SFR calculations are fairly insensitive to z, these are 
likely to have similar SFR limits to the other bursts in 
Table 2. 

Radio observations of GRB 980703 that were made 
years after the burst indicated the presence of a persistent 



radio source (Berger et al. 2001). If this radio emission is 



related to star formation, the inferred SFR (M > 5Mq) 
was estimated to be ~ 90 MQyr~^. Table 2 shows that our 
850 /im SCUBA limit is consistent with this result. Our 
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450 fim rms of 20 mJy for GRB 980703 places a slightly 
higher Hmit of SFR (M > SMq) < 550 MQyr"!. 

The fact that we are not consistently detecting bright 
sub-millimeter sources with SCUBA naturally leads to 
questions regarding the connection between GRBs and 
dust-enshrouded star formation. If GRBs are due to the 
explosions of high-mass stars, this may indicate that the 
relatively small population of extremely luminous dusty 
galaxies does not dominate the total star formation in the 
universe at early epochs. Instead, the GRBs may be pre- 
dominantly tracing slightly lower luminosity galaxies with 
SFR (M > 5Mq) - 10 - 100 Mgyr-i. 

We caution that our sample of bursts is still rather 
small. We also note that there may be an important selec- 
tion effect in the bursts studied so far. Since the current 
GRB error boxes are large, we have to wait for the de- 
termination of a counterpart at other wavelengths before 
proceeding with our observations. If a bright optical after- 
glow is found, this suggests that the absorption local to 
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